Probing Dynamics of Single Molecules: Non-linear Spectroscopy Approach 
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A two level model of a single molecule undergoing spectral diffusion dynamics and interacting 
with a sequence of two short laser pulses is investigated. Analytical solution for the probability 
of n = 0, 1,2 photon emission events for the telegraph and Gaussian processes are obtained. We 
examine under what circumstances the photon statistics emerging from such pump-probe set up 
provides new information on the stochastic process parameters, and what are the measurement 
limitations of this technique. The impulsive and selective limits, the semiclassical approximation, 
and the fast modulation limit, exhibit general behaviors of this new type of spectroscopy. We show, 
that in the fast modulation limit, where one has to use impulsive pulses in order to obtain meaningful 
results, the information on the photon statistics is contained in the molecule's dipole correlation 
function, equivalently to continuous wave experiments. In contrast, the photon statistics obtained 
within the selective limit depends on the both spectral shifts and rates and exhibits oscillations, 
which are not found in the corresponding line-shape. 

PACS numbers: 82.37-j, 82.53-k, 05.10.Gg, 42.50.Ar 
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I. INTRODUCTION 

Recently van Dijk et al [l| reported the first experi- 
mental ultra-fast pump-probe study of a single molecule 
system. Unlike previous approaches to non-linear spec- 
troscopy where only the ensemble average response to 
the external fields is resolved [1], the new method yields 
direct information on single molecule dynamics, gained 
through the analysis of photon statistics. Although the 
original experiment [l[ was conducted on a molecule 
undergoing a relatively simple relaxation process, the 
potential of combining non-linear spectroscopy with 
single molecule spectroscopy inspires many unanswered 
questions: What are the limitations of the investigation 
of fast dynamics? How does the information contained in 
these experiments differ from the information contained 
in simpler continuous wave experiments? How to design 
the external control fields, so that needed information 
on dynamics of molecules is gained? What is the finger 
print of coherence in these types of experiments, and 
how its influence on photon statistics is suppressed due 
to dephasing processes? The answers to these questions 
are important for better understanding of a wide variety 
of physical phenomena and have implication in the 
investigation of ultra-fast dynamics of molecules in the 
condensed phase, of quantum properties of light, and 
in the field of quantum information and computation 
P, H, 0, H, 0, 0, H 01 ■ Here we present a treatment based 
on the stochastic Kubo- Anderson model 0, [l^, [ll|, , 
which yields general insights on the problem. 

We consider a sequence of two identical laser pulses 
interacting with a single molecule (or an atom, or a 
quantum dot) undergoing a spectral diffusion process, 
namely a molecule whose absorption frequency is 
randomly modulated in time due to interaction with 
a thermal bath. The electronic states of the single 
emitter are modeled based on the two level system 
approximation. Most single molecules have a triplet 



state, however the life time of the triplet is much longer 
than the time scales under consideration in this paper, 
and it can be neglected. It is assumed that the pulses 
are very short compared with the inverse rate R of the 
spectral diffusion process, as well as with the inverse of 
the radiative life-time of the emitter, T. The probability 
of photon emissions during the pulse event is then 
negligible, and a pair of pulses yields two photons at 
most. Repeating the experiment many times one obtains 
the probabilities {Po),{Pi) and {P2) of emitting 0,1 
and 2 photons, where (• • •) designates the average over 
the stochastic modulation of the two level system's 
absorption frequency oj(t). In what follows we generalize 
the results obtained in our earlier publication [l3 | by: 

(i) establishing the general expressions for (Pq), (Pi) and 
(P2) in the limit of long measurement times without any 
restricting assumptions regarding the laser detuning, 

(ii) comparing the photon statistics obtained for the 
two state Kubo-Anderson and Gaussian processes. 
We show, that under certain conditions this type of 
photon statistics reveals important information on single 
molecule dynamics, information which might be difficult 
to obtain using other theoretical approaches to single 
molecule spectroscopy [11, [11 [11, [M Il3 H [ll 113, iSI . 

The general expressions for the photon statistics are 
obtained starting with the path interpretation of Mollow 
and ZoUer, Marte and Walls [13, [ll] of the optical Bloch 
equations [11]. We show, that depending on the char- 
acteristics of the stochastic dynamics and the laser field 
parameters, different types of non-linear spectroscopies 
emerge. In particular, sensitivity to the phase accumu- 
lated by the system in the delay interval between the 
pulses is found, and impulsive and selective type of spec- 
troscopies are considered in detail. The Kubo- Anderson 
spectral diffusion process [lO'l used in this work is found 
in many molecular systems [14, 19, 2(1, JlJ and may be 
easily detected, when the process is slow by means of 
the spectral trail technique. Our goal is developing gen- 
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FIG. 1: The two square laser pulses as modeled in Eq. ((Sj) on 
the background of a schematic realization of spectral diffusion 
process Lo{i). T = t\ — — 12 is the width of single pulse and 
A is the delay time between two subsequent pulses. 



state into the ground state, due to spontaneous emission 
with r designating the emission rate. Finally ult) is the 
stochastic time dependent absorption frequency of the 
system. The spectral diffusion process uj{t) is modeled 
using the Kubo- Anderson approach: 

uj{t) ^ujQ + Sw{t), (3) 

where coq is the bare absorption frequency of the single 
emitter, and d'w{t) is a random function of time llo|. 
We will assume the process 6w{t) is stationary, its mean 
is zero, its correlation function is 

{5wito + t)Swito)) ^ly^m, (4) 

where ■0(0) — 1 ^-nd ^(oo) = 0. Later we will demon- 
strate our results for a particular choice 

m = e-'^*, (5) 



eral methods suitable for detection of the wider range 
of dynamics. Some technical details of the calculations 
skipped in the main text appear in Appendixes A, B and 
C. 



II. PHOTON STATISTICS 

In our model two identical square pulses interact with 
the two level system (see Fig. 1). The state of the 
system is described by the density matrix, represented 
by the 4- vector a = (ctoc, Cgg, (Xgc, Ccg)^. Here cTcc and 
(Tgg represent the populations of the excited and ground 
states respectively and (TgcCTog describe the coherences, 
namely the off diagonal matrix elements of the density 
matrix. For the sake of mathematical convenience we use 
the following 4-dimensional basis [2^: |e) = (1,0,0,0)"^, 
which means that the system is in the pure excited state, 
|g) = (0, 1, 0, 0)'^ - the ground state, and |c) = (0, 0, 1, 0)"^ 
and |c*) = (0,0,0,1)^ describe the coherences. The dy- 
namical evolution of the density matrix in the presence 
of the external laser field is given by the optical Bloch 
equation [2^] 



a = L{t)a + Fcr. 



(1) 



The operator 



-F -i^f{t) i^f{t) 

i^fit) -inf{t) 

-injit) iflfit) iuj{t) - F/2 
inf{t) -iilfit) -iuj{t) ~r/2 

(2) 

describes the interaction of the system with the driv- 
ing electromagnetic field through flf{t), where = 
— (l/?i)dgc • Eq is the Rabi frequency with Eq - the am- 
plitude of the electric field and dge = deg the transition 
dipole moment of the two level system. The operator 
r = r|g)(e| describes the transition from the excited 



obtaining semi-analytical solution for the Gaussian 
process and analytical solution for the two-state Kubo- 
Anderson process [13 . \2^ . where uj{t) = ujo + ly or 
uj{t) — loq — V, with the rate R determining the tran- 
sition between the -I- and — states. The later is used 
to model single molecules in low temperature glasses [2l| . 

In the case of two identical square pulses the modulat- 
ing function f{t) in Eq. ^ is: 



fit) = { 



' cos {ijJLt) 



cos [u}L{t- 





< i < ti 

h <t <t2 

t2)] t2<t <t3 

h < t 



(6) 



where ljl is the laser frequency, ti = is — ^2 — T 
is the pulses duration, and A = i2 — ti is the delay 
between the pulses. In our calculations we assume: (i) 
that the system is always found in its ground state at 
the beginning of the experiment, (ii) In order to get 
meaningful measurable results one has to use sufficiently 
intense laser fields such that fiT ~ 1 - weak fields cannot 
excite the emitter even once, (iii) The pulses are short 
enough, so that the lifetime of the excited state is much 
longer than pulse's duration - FT ^ 1 (hence, SI ^ F). 
(iv) The rate R of changes of the stochastic process uj{t) 
satisfies RT <C 1 (hence, 3> i?). Assumption (iii) 
leads to the mentioned negligibility of photon emissions 
during the pulse events. The last assumption implies, 
that the time-dependent absorption frequency ult) is 
unchanged during the excitations, and will be taken to 
be uj{ti) = loq + Sw{ti) during the first pulse event and 
Lu{t2) =0^0 + Sw(t2) during the second pulse event. 

A calculation given in short in Appendix A yields the 
following expression for the probabilities of emitting zero, 
one and two photons in the limit of long measurement 
times t — > oo for a particular realization of the stochastic 
process w{t): 
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P„[lo ih) , Lu m) , 0(A)] = PCI- [u {t,),Lu ih)] + 2e-^Re { A^^h (t,) , ^ ^t,)] e^(e(A)+a.„(T+A)) J ^ (7) 
where 0(A) is the random phase accumulated during the delay interval given by 

e{A) = [ Sw{t)dt, (8) 
"'0 

and P.^^'^ [uj (ti) , uj {t2)] and A'^"^ [uj {ti) , uj {t2)] are given in Table 1, where G [uj (ti)] {i = 1, 2) is the Green function 
defined in Eqs (|9ll0p below. Note, from Eq. (O it follows that all accumulated random phase effects become negligible 
when FA 3> 1. The two summands in Eq. ([7]) describe two kinds of possible quantum trajectories leading to required 
number of photon emission events: the term P^''* [uj (ti) ,uj (^2)] we call semiclassical in the sense that it summarizes 
the paths where at the beginning of the delay the system is found in one of the pure states |e) or |g) [2^. For example, 
consider the first term of P^^^: the system starts in the electronic ground state |g), then it evolves with the propagator 
of the first pulse G without photon emissions and reaches the excited state |e), it stays in the excited state 

during the delay interval (with probability e~^'^), and afterwards the second pulse with the propagator Q [^(^2)] 
stimulates the induced emission, bringing the system back to the ground state |g) without emitting a photon. The 
second summand in the right-hand side of Eq. ([7]), A'j^°^ [oj (ti) ,uj (^2)], summarizes the contribution of the coherence 
effects (i.e. all those quantum paths where at the beginning of the delay interval the system is left in a superposition 
of the pure states). It can be shown (see Appendix A), that very strong or resonant 7r-pulses, where flT = tt, simply 
switch the state of the molecule being in |g) to |e) or vice versa. Hence, they do not excite the coherence - in such 
cases the coherent term A'^"^ [uj (ti) ,uj (^2)] vanishes. 



n 


Pf-[c.(ti),C.(t2)] 


A?,'^^[uj{h),cuit2)] 





{g\g [uj (i2)] \e){e\g [uj (ti)] |g)e-^^ + {glG [uj (t2)] \g){g\g [uj {h)] |g) 


{g\g[uj{h)]\c){c\g[uj{h)]\g) 


1 


{g\g [u; ih)] \g){e\g [lu ih)] |g) + (el^; [lu (t^)] \g){g\g [lo ih)] |g) 


{o\g[u;it2)]\c){c\g[ujih)]\g) 


2 


{e\g[u.ih)]\g){e\g[uj{h)]\g) 






Table 1: Photon statistics for two short pulses and arbitrary spectral diffusion process uj{t). 



In Table 1 the propagator of the two level system 
during the first and the second pulse events 5[a;(ti)] 
{i = 1, 2) is obtained within the rotating wave approxi- 
mation (RWA) (see Appendix A): 



G [uJ ih)] = 



where 



L 



and 



with 



RWA 
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S{ti) = 


Sl- 


6w{ti) 






-- UJL 


- UJQ 





(9) 



(10) 



(11) 



(12) 



is the detuning. The mathematical calculations were 
made with the help of Mathematica 5.0. Evidently, since 
the spontaneous emissions during the pulse events are 
neglected (il ^ F), the Green function Eq. Q describes 
well-known Rabi oscillations. 



III. INFLUENCE OF SPECTRAL DIFFUSION 
ON PHOTON STATISTICS. 



We now take the average of Pn[uj (ti) ,uj {t2) ,0{A)] 
Eq. ([7]) over the stochastic process uj{t). This procedure 
requires the knowledge of the joint probability density 
function (PDF) V[uj{ti),uj{t2),0{A)] of finding the sys- 
tem's absorption frequency uj{t) in the infinitesimal range 
near ujq + 5w{ti) at ti, near ujq -f Sw{t2) at t2, with ac- 
cumulated random phase ^(A). Then by definition the 
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FIG. 2: The matrix element {e\Q[uj{U)]\g} Eq. ((HI describing 
the probability of transition from the ground to excited state 
as a result of interaction with a vr-pulse for Q = 4r. The 
smooth curve represents (e|{?[aj(ti)] |g) for uj{t) = ujo + v with 
u = Q,. The dashed curve represents (e|C/[aj(fi)]|g) for ijj[t) — 
ujo — V with v = Q. The figure illustrates that since the two 
curves practically do not overlap for u > ^1 the probabilities 
to excite the molecule being in state = loq + v during 

the first pulse event and in the state loit^) = uq — v during 
the second pulse event are measurably different. 



average of P„ is: 



{Pn) 



d6i(A) / dw(ti) / dw(i2)x 
. Jo Jo 



X P^[u{h),L0{t2),e{^)]V[uj{h),Uj{t2),e{^)]. (13) 

Later we find exact solution for the three-variable PDF 
^(^2), 6'(A)] for the case of two-state Kubo- 
Anderson process, thus providing all essential tools for 
calculation of (P„) in the case of the telegraph noise. 
But first, we discuss several limiting cases common for 



all stationary processes. 

In Fig. 2 we plotted the probability of transition from 
the ground to the excited state (e|C/[a;(ti)]|g) Eq. ([M]) 
for two identical 7r-pulses. The smooth and the dashed 
curves represent (e|CJ[a;(ti)]|g) as a function oi 5l for 
uj{t) = uio + v and =u!o — u respectively with i/ = fl. 
The half width of (e|5[a-'(ti)] |g) is of the order of fl. 
When ;/ 3> the two curves practically do not overlap. 
Hence, assuming the absorption frequency of the system 
at the moment of the first excitation is uj{ti) = luq + v 
and of the second excitation is uj(t2) = ujq — 1^, the 
probabilities to excite the molecule during two pulse 
events strictly differ. In this selective limit the photon 
statistics is therefore very sensitive to the temporal state 
of the molecule at the moments of excitations, and the 
particular type of the underlying stochastic process has 
large importance. Later we consider this limit in detail 
for the case of telegraph noise. The opposite situation is 
the impulsive limit. 

Impulsive limit Q ^ v. - In the limit ft ^ v, 
which we call impulsive, the matrix elements of 5[ti;(ii)] 
and tj[cij(i2)] Eqs. (|54p become independent of the value 
of stochastic detuning 6w{t) at the moment of the ex- 
citation. Thus instead of the multi variable PDF 
P[uj(ti),uj{t2),0{A)] we now have to deal only with the 
one variable PDF of the phase ^(A). As a result the 
photon statistics shows an interesting relation with lin- 
ear continuous wave spectroscopy: 



lim (Pn 



lim R^'''+2e"-^Re 



where using Table 1 and Eqs. ([51)1 



(14) 



2Q2 



hm (PI 



Cla 



5l + n'cos'i^^Jl , 



2 / nx 



{61 + 



(15) 



and 



lim (pCla^ = il^^l]4gi^4 



hm (Po^'^) = 1 - lim (Pi^'^^) - hm {P^^^) 



lim A^°^ = 



I 



(16) 
(17) 

(18) 



In Eq. (fTil) the function 0(A) given by 

0(A) = (exp[i6l(A)]) = (exp[i / duj{t)dt]) (19) 



is the well investigated Kubo- Anderson correlation 
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function, whose Fourier transform is the line shape of 
the two level system according to the Wiener-Khintchinc 
theorem In conclusion, we see, that working with 

very strong laser fields f2 ^ under assumptions (iii,iv) 
we gain the same information as found in the line-shape 
in continuous wave experiments. 

Near the resonance, where (5l ~ 0, using Eqs. p4j|18p 
we find 



lim(Po) =e-rAsin4 ( 



cos 



or 



(20) 



lim (Pi) = - sin^ (nr) ( 



rA 

1 + e 5- Re 



lim (P2) = (1 - e-r^^ "^"4 



0(A)e-o(r+A) 



}■ 



sm'l— I 



(21) 
(22) 



From Eqs. (|20j|2ip we see, that for the strong 7r-pulses 
the coherent terms vanish, as mentioned. In contrast, 
for 7r/2-pulses with flT = 7r/2 the importance of the 
coherent terms, and hence, the correlation function 
(/)(A) on the photon statistics is the strongest, since 
the tt/2 pulse excites the off diagonal terms of the 
pulse-propagators Q[uj{ti)] Eq. ((54|) [2^. 



Semiclassical approximation. - The influence of the 
coherence on photon statistics in many experimental 
cases is expected to be difficult to detect. It may be 
because of the dephasing effects caused by the damping 
coefficient e~~ multiplying the coherent terms A^°^ in 
Eq. ([7]). Moreover, because of the large value of the bare 
optical transition frequency ujq the coherent paths oscil- 
late too fast to be detected (see the term e^^oiT+A) 
Eq. (O). In such cases a practical approximation is to 
keep only the semiclassical terms P^'''. Nevertheless, we 
stress, that for multilevel systems or in non optical ex- 
periments on Josephson junction coherence contribution 
is important 8]. Since the semiclassical paths are inde- 
pendent of the random phase 0{A), for the calculation of 



{Pn^'^} we need only the marginal, two dimensional PDF 
'P\uj(ti),uj(t2)]- For example, in the case of Gaussian 
noise 1281 we have: 



V[Sw{ti),dw{t2)] = 



1 



X exp 



Sw{ti)^ + 5w{t2f - 2Sw{ti)Sw(t2)^' 



21/2(1-^2) 



(23) 



Jo 



where ip = ip{\t2 — ti\) is the time dependent part of the 
correlation function in Eq. ([4]). Once the two-time PDF 
V[Lu{ti),uj{t2)] is known, (P„) within the semiclassical 
approximation is: 



(pCia^^Z / P^'-[u;{h),u;{t2)]V[u;{h),u;{t2)]du;{h)du;{t2). 

(24) 

In Appendix B we give the explicit semiclassical ap- 
proximation for the two state Kubo- Anderson model 
Eqs. (|58ll57p . Later in Fig. 3 we compare these results 
with similar results for the Gaussian noise, which was 
solved semi-analytically with the help of Mathematica 
5.0. In the both calculations the same correlation 
function Eqs. (|4I5P was used. 



IV. TWO STATE PROCESS: EXACT SOLUTION 

Now we obtain the exact solution for the two state 
Kubo- Anderson Poissonian process, where the absorp- 
tion frequency of the system jumps between the -|- and 
— states, i.e. uj{t) = loq ± v. We denote the ini- 
tial state, during the first pulse with uj{ti) = + or 
Lo{ti) = —. Similarly, the final state at the second pulse 
is w(t2) = + or — . Since the random phase is now given 
by /S.9 = u{T^ — T~), where are occupation times 
in states + and — [23| obeying A — + T~ , the joint 
PDF 'P[uj{ti),uj{t2),9{A)] can be found from the joint 
PDF h[uj{ti),uj{t2),T^] of finding the system in state 
Lu{ti) — ± during the first pulse, state uj{t2) = ± during 
the second and with the occupation time between the 
two pulses. In this case, where the random process takes 
only discrete values, (P„) Eq. (fT3)) takes the form 



(Pn) ^['^(ii),^(i2)]Pf"[c^(ti),c.(^2)]+2Reje^'^''(^+^)-^^^oh[^(^^)^^(^^)] / e'-(^^^-'^^h[u;{h),u;it2),T+] dT+ 

(25) 

where 

1 I -2RA 1 _ -2_RA 

7'[±,±]-— ^ , ^[±,T] = 1 (26) 
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are the probabilities of finding the particle initially in state Lu{ti) = ± and finally in state ^(^2) = ±, which are easy 
to obtain from Poissonian statistics. In the integrand of Eq. (P5|) one can recognize Laplace transform r+ —2iv of 
h [aj(ti), [^(^2), Simple rearrangement leads to: 

U{tl),uj{t2) = ± 

(27) 

where h[io{ti),Lu{t2),—2ih'] is the Laplace T+ —^ —2iv transform of h[u>{ti),u!{t2),T^]. The procedure of derivation 
of h[u;{ti),uj{t2),—2iiy], based on the renewal processes theory [13, [2§|, is given in detail in Appendix C. Here we 
present the final results: 



h[±,±,-2iL'] = 



coshlAV^R^T^I ± 



h [t, ±, -2w] = e 



'sinh(AVi?^ - 1/2) 



(28) 



(29) 



2Vi?2 - 1/2 

Using Eqs. (|26II29I54)) and Table 1 the calculation of (P„) is straightforward. 

Semiclassical selective limit. - Now we consider selective limit for the two state Kubo- Anderson process within the 
semiclassical approximation. From the exact solution Eqs. (j56ll58l) . considering two opposite situations Q ^ Sl + 1^ 
and Q <ti Sl — V or fl 6l — and f2 ^ + we find: 



(l+e-2fl^)r!2 sin2 



( lim Pi*^'^) = 



1W 



2[{\5L\-iyr + n^] 



(|<^l| -1^)'+^^' COs2 






)2 +172 



2[{\SL\-iyr + n^] 

(l + e-2«^)(l-e-rA)r!4 sin^ 



4([(|'5l|-^^)2+172]^ 



lim Pq^'^) = 1 - ( hm Pi^'^> - ( lim P^^"") 



V»r2 



(30) 

(31) 
(32) 



Notice, that P^'^ia depends on A only trough e 2-RA jjgjj^g^ using Eq. ((30|) it is easy to measure the dynamics of the 
molecule. 



r 



When 5l = ±1^ from Eqs. (|30ll32p we find: 



hm (P,C'^) = sin2(f7r)+^^ sin2 ('^ 

|5i,|=i/»o 8 2 \ 2 



(33) 



lim {P^^^) 



(l + e-2^^)(l~e-^^) . 4 /ra; 

4 2 



(34) 



Note, that Eqs. p3l34p exhibit Rabi oscillations. Apply- 
ing Eqs. (|33l34p to the 7r-pulses gives: 



lim (Pf = i (1 - e 



1 



-2_RA\ 



(35) 



and 



hm (P„c>^) = ± (1 + e-^^^) (1 + e-^^) . (37) 
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FIG. 3: Photon statistics for two identical 7r-pulses, where A — SF"^, _R = F = 1. In the first row: the probabihties of emitting 
a. b. 1 and c. 2 photons in the case of the selective limit v — 50F, f2 — 20F. The exact solution for the two-state 
Kubo- Anderson process Eqs. (I56j|58|l is represented by the smooth curves. The dot-dashed line represents the semiclassical 
selective limit approximation Eqs. (I30I32|I for the two state process. And the dashed curves show the Gaussian process.; In 
the second row: the probabilities of emitting d. e. 1 and f. 2 photons in the intermediate situation v = 20F, — 20F.; In 
the third row: the probabilities of emitting g. h. 1 and i. 2 photons in the case of the impulsive limit v = 2F, f2 = 20F. The 
dot-dashed line represents the semiclassical impulsive limit approximation Eqs. (|15I17|) . (The curves representing the exact 
solution for the telegraph and Gaussian noises, and the impulsive limit approximation coincide). Pay attention, the graphs are 
scaled differently. 



The results of Eqs. p5ll37p make perfect physical sense. 
For example a single photon may be emitted only if 
the absorption frequency uj{t) is found once in the + 
state and once in the — state, since only one of these 
states is in resonance with the laser. The system, 
which is in the ground state at the beginning of the 
experiment, gets excited exactly once, and nothing 
interrupts the spontaneous emission process. Hence, 
(Pi) = +'P[~,+]. Note, that for the strong 

TT-pulses in the case of the impulsive limit we have 
(Pi) (sec Eq. (Ell)). 



V. DEMONSTRATION OF RESULTS 

In Fig. 3 we plotted the semiclassical parts of (Po), 
(Pi) and (P2) for the two state and Gaussian processes 



for two identical 7r-pulses. We see the transition from 
the selective limit (first row) to the impulsive limit 
(third row), where the graphs corresponding to the two 
processes visually coincide. The graphs in Fig. 3 clearly 
provide the information on the spectral shifts ly. Finally, 
it is worth noticing, that the behavior of photon statis- 
tics corresponding to the two-state process is oscillatory 
in 6l, while in the case of the Gaussian noise it is not. 
The origin of this effect follows from the fact, that the 
matrix elements of the propagators G [uJ (ti)] Eqs. (I9l54p . 
and hence the probability of n photon emission events, 
P„ are sinusoidal functions of Sl — Sw(t). Therefore 
in the case of discrete dichotomic noise the phase of 
the sinusoidal functions we are summing up takes only 
two values ±i/. However, in the case of the continuous 
Gaussian distribution we integrate over a continuous 
range of phases. Thus, as we approach the selective limit 
in the case of Gaussian noise, the oscillations in 6^ are 
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FIG. 4: The plot of {Pi) for the two state Kubo- Anderson 
process as a function of the delay A and spectral shift in 
the case of pure resonance with the + state Sl = v for two 
identical 7r-pulses with Q = 60r, R = 3T. The lower graph is 
the contourplot of the upper three-dimensional graph. 

destroyed by averaging. 

In Fig. 4 we plotted the three dimensional (above) and 
contour (below) graphs of (Pp^^) for the two state Kubo- 
Anderson process as a function of the delay interval A 
and the spectral shift v. The graphs were obtained for 
the case of two 7r-pulses for a constant values of = 60r, 
R — 3r, r = 1 (hence, RT ~ 0.157) in pure resonance 
with the + state Sl = v- These graphs clearly show, 
that the photon statistics of the pump- probe set up yields 
the dynamical information on the rate R. However, for 
1/ <C this information cannot be obtained, since (Pp'**) 
stays visibly unchanged along the A expected. 
For larger shifts the changes along the A-axes become 
more and more significant, and (Pf^''^) can be used to 
measure R. 



(c) 1/ n R — selective pulse, slow process, 
while all the other combinations violate our assumptions. 

As mentioned, the direct dependence of the photon 
statistics on the shift v and the rate R becomes unde- 
tectable when the strength of the laser field reaches the 
impulsive limit. We saw (see Fig. 2), that for < the 
information on the stochastic process encoded in photon 
statistics becomes equivalent to the information provided 
by line-shape in continuous wave spectroscopy. Since 
il ^ R we see, that the effectiveness of the pump-probe 
technique is restricted to slow processes v > R. Thus, 
the assumptions RT <C 1 and ^ 1 determine the 
limitations of the measurement of stochastic dynamics, 
and to get information on wider range of dynamics our 
method must be improved by removing these conditions. 
In our future work we plan to focus on the investigation 
of processes, where RT ^ 1 is not fulfilled. 



VI. FAST MODULATION 

Finally, let us consider an interesting case of the fast 
modulation limit i? ^ c», i/ — > c» such that R ^ v, 
where the motional narrowing effects take place. From 
inequalities (a, b, c) follows, that the fast modulation 
limit must be impulsive limit as well. Therefore, Eqs. (j20I 
hold with 

lim 0(A) = exp(-rsDA/2), (38) 

with 

TsD = i^VR, (39) 

which means that (Fsd + r)/2 is the renormalized decay 
rate which damps the coherent terms. Note, that Tsd is 
a measurable physical observable determining the width 
of the line shape in continuous wave spectroscopy |10j . 

The results of Eq. (|38l39p are obtained straightfor- 
wardly for the two state Kubo- Anderson process using 
our exact solution Eqs. (|28l29p . We calculate the limit 
R ^ 1/ in such a way, that Fsd — / R is finite. For 
example considering Eq. (P^)) we have: 



Measurement Limitations e~^^^~ 



Combining the requirements UT ^ 1 and RT <C 1 



I V i?2 — g V 



(40) 



implies 17 > i?. Hence, only the three following hmiting Using y^l - - 1 - ^ in the exponent of Eq. gOl) we 

obtain: 



cases may occur: 



(a) 17 >• » i? — impulsive pulse, slow process „ g-A(_R-ii. 

(6) 17 >• i? >• — impulsive pulse, fast process r^v^^~^'~^' ^^^^ ^ l>i 4 \ 
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~e ^ i^M = -e exp(- — — ), (41) 

provided _RA ^ 1. Exactly the same procedure applied 
to Eq. leads to similar result. Inserting the result 
Eq. (|¥T|l into Eq. (|77|) cancels the oscillating term e'^". 
Therefore, the stochastic phase in the limit of fast 

modulation reduces to exp {^—%^^- It is easy to show, 
that Eq. is valid as well for the Gaussian process 
under investigation, with Tsd — j^R. The result 
Eqs. ((38|) obtained in this section for the telegraph noise 
exhibits the fast modulation limit behavior general for 
all Markovian processes • 



improved by removing these conditions. In particular, 
we plan to investigate photon statistics for processes 
with KT ^ 1. In addition, our techniques can be easily 
modified for the investigation of relaxation processes [1| , 
where the frequency shift may be large, thus allowing the 
measurement of dynamics up to the order of pico-seconds. 

5. Finally, in the fast modulation limit where the 
stochastic transitions rate R is much larger then spec- 
tral shifts I/, the influence of the underlying spectral 
diffusion dynamics on the photon statistics reduces to 
the well-known Kubo-Anderson correlation function 
coefficient exp {^—^T^sd)^ which damps the contribution 
of the coherent terms. 



VII. SUMMARY 
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VIII. APPENDIX A 

Photon Statistics. Starting with (2^ . [23l | an interpre- 
tation of the optical Bloch formalism yields a tool for the 
calculation of photon statistics. The formal solution to 
Eq. Il]) may be given by the infinite iterative expansion 
in f M: 



1. Unlike the line-shape in continuous wave experi- 
ments, the photon statistics obtained from pump-probe 
set up depends not only on the spectral shifts v, but 
also on the rate R of the spectral diffusion process and 
exhibits oscillations not found in a line-shape. 

2. In the limit of the impulsive pulses the spectral se- 
lectivity is lost. The information on the spectral diffusion 
parameters, contained only in the phase accumulated 
during the delay interval, is then equivalent to the 
information gained from a line-shape in continuous wave 
spectroscopy. In optics using pump-probe technique this 
information is expected to be difhcult to detect, because 
of the huge bare optical frequency too and dephasing 
effects. It will be interesting to investigate the coherence 
effect in multilevel systems. 

3. In contrast, in the selective limit, where the laser 
field is weak compared to the deviations of the spectral 
shifts, the photon statistics offers full information on 
spectral diffusion parameters. 

4. We have shown, that the limitation of the mea- 
surement of stochastic dynamics is determined by 
the conditions RT ^ 1 and flT ~ 1, restricting the 
effectiveness of pump-probe set up only to slow processes 
I' ^ R. Hence, to get better results our method must be 



= O)a(o) + / dtig{t,ti)tg{ti,0)a^o) + 



dt2 / dtig{t,t2)tg{t2,ti)tg{ti,o)a 



(0) 



(42) 



where tT(o) is the initial condition, and the Green function 
describing the evolution of the system in the absence of 
spontaneous transitions into the ground state (i.e. with- 
out r ) is 



g{t,t')=fexp 



L{ti)dti 



(43) 



where T is the time ordering operator. 



Each term in the expansion Eq. describes the 

propagation of the system emitting exact number of pho- 
tons: for example the first term does not include F at all, 
and describes a process where no photons are emitted, 
the second term corresponds to the processes where only 
one photon is emitted and so on. Therefore, cr|"j'' defined 



as: 



(") „ 

(t,0)'^(0) 



(44) 
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describes the conditional state of the system at the mo- 
ment t, provided that n photon emission events occurred 
in the time interval (0,t), and 



U, 



(») 
(*,*') 



dtn ■ ■ 



dti g{t,tn)r---rg{tut') (45) 



is called the n-photon-propagator. 

The main equation for calculating the probability of n 
emission events up to time t is: 



P„W = ((e| + (g|)a;")) = ((e| + {g\)Ull['\a^o)) , (46) 



(«) 



'(t) 



(t,0)l 



which is simply the trace of the density matrix condi- 
tioned by n emission events. We showed in [2^ that 
in the case of two separated pulses the total n-photon- 
propagator acting from t = wp to t may be written 
as: 



rr(") _ 7-r("-a-/3-7)7-r(7) rrl^M; tt'. 



r(/3) 



r(") 



(47) 



where the superscripts a, (3 and 7 are all non-negative in- 
tegers leading to n emission events (i.e. n—a—(3—^ > 0). 
The Einstein's summation rule from to n must be 
applied to every superscript appearing twice. Eq. (I47p 
means that the n-photon-propagator acting in (0, t) can 
be decomposed into the sum of all possible products 
of the a-photon-propagator acting in (0,ti), /3-photon- 
propagator acting in (ii,t2), 7-photon-propagator acting 
in (^2, ^3) and {n ~ a — (3 — 7)-photon-propagator acting 
in (^3,^)- 

The propagators acting during the delay interval and 
after the second pulse, where the laser is off and $7 = 0, 
may be found immediately: 



/e-rA 



(0) 

(ti+A,ti) 



1 



e''^"'^" 

















o\ 

e-rA 


_Q n n n/ 



(48) 



(49) 



and 



U, 



(ti+A,ti) 



for n > 1. 



The propagators acting in (i3,i), where t 
given by the limit A — ^ cx3. 



(50) 



00, are 



To obtain photon statistics for short pulses in the 
summation Eq. (|47p we pick up only those processes, 
where the propagators acting during the pulse events 



are the zero-photon-propagators, i.e. U, 



and U, 



(0) 



(0) 

(ti,0) 



^ ^[^^(^2)]- We substitute Eqs. (gZl Si 
|49| into Eq. (|46)l . and inserting the closure relation 
X]j=e g c c* ~ 1 between each two propagators 

obtain Eq. ([7]) of the article. 

The calculation of the matrix elements ^/[^(ti)] is 
made using RWA. Applying RWA to the optical Bloch 
equation ([T|) consists of neglecting fast oscillating non- 
resonant terms [l^l . As a result the following equation is 
obtained: 



•RWA 



= L 



RWA 



{t)a 



RWA 



RWA 



(51) 



where 



L 



RWA 














-in 


iO 


1 


2 


, iO. 
\ 2 


-iO 


2 



—in in \ 

2 2 



-jp. 
2 

2 



2 
in 
2 



-^-i(5(i) 





(52) 

and 5(t) Eq. (fTTj) is the detuning at moment t. As dis- 
cussed in the article, the detuning is assumed to be con- 
stant during the pulse. Hence, in the new representation 
the calculation of the Green function is simple: 



g\u:{U)\ = exp [rL«^A[c^(i,)]] 



(53) 



Since the interaction time goes to zero FT — > 0, the spon- 
taneous emission effects during the pulses are completely 
suppressed, and we find: 



sin^ 



(e|e[c.(iO]|g) = (&\Q\^{Um = 1 - (&\Q\^{m%) = 1 - {elGHU 



= -{c\g[Lu{ti)]\e) = ~{c*\g[u;{U)]\er = -{e\gHti)]\c) = -{e\g[Lj{U)]\c*y = 
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f7r 



2[(5z,-5u;(ij))2 + r!2] 

For pulses satisfying — — <C 1, i.e. resonant or impuslive pulses, from Eq. (|54p we find: 



(54) 



/ cos2 ^ sin^ ^2: 



\ 



■sinOT 
'* 2 
; sin QT 



COS 



2 or 

2 

sin f2T 



; sin QT 



■ sin OT 
2 

,■ sin SIT 
' 2 

os^f: 

•2 nr 



■ sin SIT 
' 2 

„■ sin nr 

'* 2 

lin^if 

2 SIT 
2 



^ COS^^ / 



(55) 



Substituting QT = tt into Eq. ([55]) we see that the off-diagonal terms giving rise to coherent terms vanish. More 
detailed discussion on RWA and its application to the calculation of the matrix element may be found in [25j . 



IX. APPENDIX B 



Semiclassical approximation for the two state Kubo- Anderson process. Here we present the exact expressions for 
{Pn^^) (n=0,l,2) for the case of the two state Kubo- Anderson process. Using Eqs. ([26]) . Table 1, the matrix elements 
Eqs. dSll) and Eq. ^ we find: 



Cla\ 



(1 



-2RA\ 



sin^ 



nr 

2 



1 



1W 



cos2 



or 

2 



1 



sin^ 



(Sl + yf + cos2 



{5l + y? + QL" 



{5l + yf + 



(1 - e-2^*^) 



fl2 sin2 



^Jl 



~sT2— 



{Sl + vf + Q? cos2 



{5l + i^Y + 



sin^ 



2 



— fir~ 



{5l -vf + VL' cos2 



QT 
2 



— TP — 



(Jl + 1/)2 + 



(5i - l/)2 + r!2 



(56) 



-2fiA 



)(1 



Q'^ sin^ 



f^A/l-t-^ 



f72)^ 



f]" sin^ 





f}2)2 



(1 



-2_R.A 



)(1 



-rA\ 



sin^ 




sin^ 


¥^1 + 




(-5 


L - ^)2 + ^2 


- X 


5l + ^^)2 


+ 02 



(57) 



and 



(58) 
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The only approximation made in the semiclassical Eqs. (|56ll58p is neglecting the spontaneous emission effects during 
the pulse events and applying RWA within the calculation of the matrix elements. Since these expressions involve the 
both V and i?, they may be used for the determination of t he spectral shifts and rates. 



X. APPENDIX C 



independently of the number of jumps. 



Derivation oj h\Lo{^\)^ 10(1-2) ^—1iv\. The two-state 
Kubo- Anderson process has been investigated exten- 
sively by many authors, in particular different techniques 
for calculation of the occupation times were proposed 
[13, m. Here we follow the method used in [li]. Provid- 
ing the molecule was in the -I- state in the beginning of 
the first pulse and in the — state in the beginning of the 
second, so that exactly m = 2k + 1 where k = 0, 1, 2 • • • 
jumps occurred, it is possible to show [2^, that the dou- 
ble Laplace transform of /m,A(r+), the PDF of the occu- 
pation time in the upper state T+ for a fixed A, is given 
by: 



i-x(g) 

2s - 



(59) 



where 



X{s+u) 



(60) 

are Laplace transforms of %(''') - the PDF of the waiting 
times between the subsequent jumps. The factor 2 in the 
nominator of Eq. (|59p is due to the probability to find the 
molecule in the -I- state in the beginning of the first pulse, 
which is equal ^ because there is symmetry between the 
states and we assumed stationary process. Summing up 
k from to 00 we have: 



s 1~x(s + m)x(s)' 



(61) 



where hs{+, —,u) is the double Laplace transform of the 
PDF of the occupation time in the -I- state, providing 
the system was in + state in the beginning of the first 
pulse and in the — state in the beginning of the second 



Now applying Eq. ([CTjl to the Poissonian process: 

X(r) = Re-"\ (62) 

which corresponds to the correlation function iP{t) given 
by Eq. ([5]), and substituting u = —2iv we find: 



-2iv] 



R 



2s{2R+s)-4i{R + s)v 



(63) 



Finally the inverse Laplace transform s — s- A leads to 
Eq. dMl) of the article: 



h[+,-,-2ii^] 



_^ .sinh [AVm 



R 



2Vi?2 - 1.2 

Due to the symmetry of the stochastic process 
h[+,—,—2ii'] — h[—,+,~2ii']. The procedure of 
derivation of Eq. (|28|) for h [±, ±, ~2ij/] is similar. 

Note, the marginal probabilities V [ijj{ti),uj{t2)] of find- 
ing the molecule in state uj{ti) — ± during the first pulse 
and in state ^(^2) — ± during the second pulse, inde- 
pendently of the values of occupation times, are easy to 
derive from Eqs. (|28l29p by setting v = Q which is equal 
to integrating out the T+ from h [w(ii), ^(^2), T^]- Thus, 
the following obvious expressions are obtained: 



P[±±] = /i[±,±,0] = 



1 + cxp(-2i?A) 



and 



1 - exp(-2i?A) 
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